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ABSTRACT 


Orifice  and  louver  sound  power  spectra  are  investigated,  using  an  intensity 
probe,  at  various  wind  speeds  in  a  low  noise,  semi-anechoic,  subsonic  wind  tunnel 
for  free  stream  velocities  below  50  meters  per  second.  The  radiated  noise  is 
created  by  turbulent  flow  over  various  orifice  and  louver  geometries  which  are 
flushed  mounted  into  the  wall  of  a  long  duct.  Five  orifice  samples  of  rectangular 
shape  and  various  transverse  dimensions  as  well  as  four  louver  samples  with 
multiple  rectangular  and  circular  orifices  are  tested.  Also  investigated  is  the 
effect  of  the  leading  and  trailing  edge  angle  on  the  radiated  sound  power.  The 
scaling  laws  of  the  excitation  frequencies  and  the  speed/power  laws  are  presented 
for  ratios  of  the  boundary  layer  thickness  to  the  transverse  orifice  dimension 
from  1.01  to  4.29.  A  detailed  theoretical  model  is  developed  for  rectangular 
shaped  aperture  orifices  and  louvers  based  on  the  work  by  Ffowcs  Williams  [9] , 
Nelson  [11],  and  Corcos  [12].  The  model  describes  the  predominant  variables 
which  effect  the  radiated  flow  noise.  To  ascertain  the  validity  of  the  model, 
it  is  used  to  collapse  the  experimental  power  spectra  for  two  orifice  and  one 
louver  geometry.  The  model  showed  excellent  agreement  with  the  actual  sound 
power  measurements.  Directivity  plots  are  also  presented  to  further  describe 
the  orifice/louver  source. 


1.  INTRODUCTION 


1.1  Background 

Turbulent  flow  over  an  orifice  or  louvered  openings  can  generate  high  levels  of 
broadband  noise.  This  phenomenon  has  been  known  for  some  time  [lj  and  early 
investigators  like  Dunham  |2],  Rossiter  [3]  and  East  [4]  attempted  to  identify 
and  control  the  high  vibrational  levels  and  possibly  the  radiated  sound  levels 
caused  by  the  interaction  of  the  flow  over  these  openings.  Recent  work  by  Elder, 
Farabee  and  DeMetz  [5]  concentrated  on  the  area  of  cavity  resonance  where  the 
opening  provides  the  feedback  path  between  the  closed  acoustic  cavity  and  the 
source,  the  turbulent  flow.  When  the  typical  dimension  of  the  cavity  is  large 
compared  to  the  acoustic  wavelength,  the  source  can  excite  various  room  modes 
of  the  cavity.  If  the  turbulent  vortex  shedding  frequency  coincides  with  one  of 
these  room  modes,  then  a  strong  feedback  will  occur,  resulting  in  an  acoustic 
resonance  of  that  mode. 

If  the  cavity’s  volume  is  further  increased,  multiple  room  modes  will  be  ex¬ 
cited  and  the  strong  feedback  mechanism  will  be  lost.  This  results  in  a  much 
broader  acoustic  response  and  at  lower  amplitude  levels  than  the  resonant  case 
which  is  typically  characterized  by  large  amplitudes  at  very  discrete  tonal  fre¬ 
quencies.  In  practice,  the  detection  of  resonating  cavities  is  often  easy  due  to 
their  strong  discrete  tones.  However,  once  these  resonances  have  been  corrected, 
or  presumably  prevented  in  the  design,  identifying  the  cause  of  the  remaining 
flow  noise  becomes  much  more  difficult.  Possible  contributors  to  the  broadband 
flow  noise  are:  (1)  the  turbulent  boundary  layer  direct  transfer  of  energy  to  the 
structure  of  the  enclosed  acoustic  space,  (2)  the  power  flow  through  the  opening 
(open  window  effect)  caused  by  the  acoustic  field  created  by  the  upstream  and 
downstream  turbulent  sources  in  the  boundary  layer,  (3)  the  power  flow  into 
the  acoustic  space  caused  by  the  source  created  in  the  opening  due  to  the  driv¬ 
ing  aerodynamic  field  and  the  shear-layer  edge  interaction  with  the  downstream 
edge  of  the  orifice  and  (4)  the  vibrational  excitation  of  the  internal  structure 
caused  by  the  orifice  source.  The  amount  of  noise  which  the  latter  contributes 
depends  on  the  details  of  the  internal  structure  and  the  orifice/louver  attachment 
design.  Therefore  to  predict  the  contribution  from  internal  structure  vibration, 
one  must  first  determine  the  orifice  source. 
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1.2  Content 


The  contribution  of  the  first  three  possible  contributors  to  the  total  power 
radiated  into  a  closed  acoustic  space  is  analyzed  based  on  the  results  from  sev¬ 
eral  experiments.  It  is  shown  that  indeed  the  orifice  source  provides  both  the 
spectral  character  and  the  predominant  amount  of  the  total  sound  power.  By 
analyzing  and  comparing  the  results  with  a  derived  equivalent  source  model,  the 
nature  of  the  source  and  the  parameters  which  affect  its  strength  and  spectral 
character  are  identified.  The  analysis  comes  directly  from  experimental  mea¬ 
surements  taken  in  a  subsonic  wind  tunnel  at  the  M.I.T.  Acoustic  and  Vibration 
Laboratory. 
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2.  SOUND  INTENSITY  MEASUREMENTS 


2.1  Why  Measure  Sound  Intensity? 

To  describe  an  acoustic  source,  one  must  measure  a  quantity  which  is  not 
dependent  on  the  environment  or  the  surroundings.  The  quantity  most  often 
identified  to  describe  a  source  is  its  power.  Knowing  the  source’s  power,  its 
interaction  in  any  known  environment  or  any  known  surrounding  can  be  deter¬ 
mined. 

Until  recently,  most  all  experimental  measurements  to  describe  an  acoustic 
source  were  accomplished  by  measuring  sound  pressure.  This  quantity,  however, 
was  the  result  of  a  source’s  interaction  with  a  specific  environment  and  utilizing 
the  experimental  results  in  other  applications  was  often  difficult.  One  technique 
that  was  developed  to  overcome  the  difficulty  in  relating  one  environment  to 
other  was  to  isolate  the  source  of  interest  and  to  simulate  a  free  field  environment 
in  an  anechoic  chamber.  Here  the  radiated  sound  pressure  in  air  was  measured, 
then  related  directly  to  the  sound  intensity.  This  technique  has  as  its  biggest 
limitation,  the  availability  of  an  effective  anechoic  chamber,  particularly  for  u^e 
in  low  frequency  analysis  at  the  low  frequencies. 

Another  technique  is  to  measure  the  intensity  directly  utilizing  a  sound  in¬ 
tensity  probe.  The  theory  of  intensity  measurements  comes  directly  from  two 
fundamental  properties  of  acoustics.  Those  are  the  conservation  of  momentum 
and  the  definition  of  intensity.  Conservation  on  momentum  states  that: 

pdvr/dt  —  ~dp/dr  (2.1) 

which  can  be  solved  for  the  radial  component  of  the  particle  velocity,  vr 

vr  -  -(1/p)  j  dp/dr  dt  (2.1a) 
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A  finite  difference  approximation  is  used  to  calculate  the  pressure  gradient 
using  two  identical  microphones  by 


dp/ dr  =  lim  (Ap/Ar)  «  (pA  -  ps)/Ar 

Ar-*0 


(2.2) 
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where  A  and  B  correspond  to  channels  A  and  B  for  the  respective  microphone 
inputs  and  A r  is  the  finite  microphone  spacing.  The  instantaneous  pressure  is 
taken  as  the  mean  value  of  the  two  pressure  signals, 

P  «  (Pa  +  Pb)/ 2.  (2.2a) 


The  finite  difference  approximation  is  valid  provided  the  separation  distance, 
Ar,  is  small  compared  to  the  acoustic  wavelength.  Thus,  a  finite  difference  ap¬ 
proximation  error  becomes  significant  at  high  frequencies.  For  a  plane  sinusoidal 
wave,  propagating  along  the  axis  joining  the  two  microphones,  the  finite  differ¬ 
ence  approximation  assumes  the  free  field  phase  between  the  two  microphone 
positions  to  be  equal  to  fcAr,  whereas  the  actual  measured  free  field  phase  is 
sin(kAr).  Therefore,  the  measured  intensity,  /',  is  related  to  the  actual  intensity, 
/,  by 

/'//  =  sin(kAr)/kAr  (2.3) 


and  the  resulting  measured  intensity  underestimates  the  actual  intensity  by  Le 
at  high  frequency  by 

Lt  =  10  log(sin(fcAr)/fcAr).  (2.3a) 


All  the  intensity  measurements  in  this  report  were  obtained  with  a  spacing  of 
Ar  =  50  mm.  Thus  the  finite  difference  error  for  this  spacing,  at  the  upper 
analysis  bandwidth  frequency  of  /  =  1600  Hz,  was  Le  =  —1.7  dB. 

From  the  definition  of  intensity  and  by  substituting  equations  2.2  and  2.2a, 
the  measured  sound  intensity  becomes 

Ir  =  pvT  «  -[(pA  +  pa)/2pArj  j ( pA  -  pB)dt.  (2.36) 


By  measuring  the  intensity  over  a  control  surface  which  surrounds  the  source, 
the  power  can  be  calculated  from: 


(2.4) 


P  =  j^Ird A 


The  advantages  of  measuring  intensity  directly  are:  (l)  the  effects  from 
diffuse  noise  contaminating  the  experimental  results  are  eliminated  since  the 
integral  of  this  noise  intensity  around  any  closed  control  surface  is  zero,  (2) 
the  sound  power  calculated  describes  the  acoustic  source  in  a  form  which  is 
readily  usable  to  design  engineers,  (3)  measurements  can  be  made  in  the  source’s 
near  field  and  (4)  there  are  no  restrictions  on  the  shape  and  size  of  the  control 
surface.  A  simple  yet  effective  analogy  to  describe  the  rationale  in  measuring 
sound  intensity,  thus  power  directly  follows: 

•  In  attempting  to  design  a  heating  system,  an  engineer  solicits  specifica¬ 
tions  from  two  companies.  One  company  describes  its  heating  system  by 
stating  the  electrical  power  and  the  thermal  efficiency  of  the  unit.  The 
other  company  describes  the  temperature  at  various  locations  inside  an 
experimental  enclosure  at  steady  state  conditions. 

It  is  obvious  in  this  example  that  describing  the  temperature  field  within  one 
environment  is  practically  useless  to  the  design  engineer  whereas  knowing  the 
power  and  thermal  efficiency,  the  designer  can  determine  the  interaction  of  this 
source  with  his  particular  environment. 

In  the  experiments  which  follow,  every  attempt  was  made  to  utilize  the  ad¬ 
vantages  of  measuring  intensity  directly  and  simulating  a  free  field  environment 
by  taking  the  measurements  in  a  semi-anechoic  chamber,  thus  minimizing  anv 
enclosure  interactions  which  could  affect  the  radiated  sound  power. 

2.2  Ideal  Versus  Practical  Sound  Intensity  Analyzer 
2.2.1  Active  and  reactive  sound  fields 

Sound  fields  exist  in  two  parts:  an  active  part  which  transports  sound  energy 
by  virtue  of  the  fact  that  the  sound  pressure  and  the  particle  velocity  are  in  phase 
and  a  reactive  part  which  stores  sound  energy  and  in  which  the  sound  pressure 
and  particle  velocity  are  in  quadrature.  A  sound  field  is  said  to  be  more  reactive 
the  more  energy  is  stored  relative  to  that  transported.  Sound  field  reactivity 
is  defined  as  the  difference  between  the  measured  sound  pressure  level,  Lp,  and 
the  maximum  intensity  level  Li. 
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2.2.2  Ideal  sound  intensity  analyzer 


An  ideal  sound  intensity  analyzer  would  measure  only  the  active  part  of 
the  total  sound  field  and  would  indicate  sound  intensity  values  approaching 
0  watts/m2  in  a  highly  reactive  (highly  diffuse)  sound  field. 

2.2.3  Practical  sound  intensity  analyzer 

A  practical  sound  intensity  analyzer  can  never  display  such  a  low  intensity 
level  due  to  the  small  but  ever  present  phase-mismatch  between  the  two  micro¬ 
phone  channels.  It  is  this  phase  mismatch  which  accounts  for  the  lower  limit  in 
the  dynamic  range  of  a  practical  sound  intensity  analyzer. 

2.3  Determining  the  Limits  of  Error  for  the  Measured  Intensity 

Errors  are  introduced  in  any  practical  intensity  measurement  by: 

1.  the  influence  of  the  phase  mis-match  in  the  probe,  <t>m  and  the  analyzer, 

<t>a, 

2.  the  influence  of  the  phase  of  the  sound  field  of  known  reactivity,  4>f,  and 

3.  using  a  particular  spacer  between  the  two  microphones. 

The  limits  of  these  errors  can  be  determined  using  a  phase-reactivity  nomo¬ 
gram  [see  reference  6j.  This  nomogram  graphically  determines,  <f>/  for  a  given 
measured  reactivity  (L,  —  Lp)  and  a  particular  spacer,  A r,  according  to  the 
relationship: 


U  -  [Lp  +  0.16)  =  10  log[l000/Ar  x  C//  x  4>,/360]  (2.5) 


The  phase  of  the  sound  field,  4>f ,  together  with  the  phase  of  the  probe  and 
the  analyzer,  —  —  and  0/  +  0a  +  <£m»  set  the  limits  of  error  in  the  measured 

phase  at  a  particular  frequency.  From  these  phase  limits,  the  nomogram  provides 
the  limits  of  error  in  dB  for  the  measured  intensity. 
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2.4  Sound  Intensity  Probe 


The  specific  details  on  the  operation  of  the  intensity  probe  used  in  this  ex¬ 
periment  are  contained  in  reference  (7).  Briefly,  the  intensity  probe  consists 
to  two  B  &  K  phase  matched,  1/2"  condenser  microphones,  type  4177/4165, 
separated  by  a  spacer  50  mm  in  length  and  on  a  common  axis.  One  micro¬ 
phone  is  a  reference  while  the  other  microphone  provides  the  pressure  gradient 
measurement. 

2.5  Intensity  Measurements 


2.5.1  Calculating  intensity 

Sound  intensity  measurements  were  analyzed  using  a  B  &  K  2032  Spectrum 
Analyzer.  The  sound  intensity  was  calculated  according  to  the  relationship: 


I{k)  =  -(l/u&rp)Im\GAB{k)}  (2.6) 


where  Ar  is  the  microphone  spacing  in  meters,  A  corresponds  to  Channel  A 
which  is  the  reference  microphone  input,  B  corresponds  to  Channel  B  which  is 
the  second  microphone  input,  k  is  the  frequency  index  and  p  is  the  density  of 
air,  calculated  according  to  the  relationship: 

p  =  0.348534p/(273  +  T)  (2.7) 


where  p  is  the  pressure  in  mbar  and  T  is  the  temperature  in  degrees  centigrade. 

The  cross  spectrum,  GAB{ fc),  is  the  complex  ensemble  average  of  the  complex 
products  of  the  complex  conjugated  one-sided  instantaneous  spectrum  GA(k) 
and  the  one-sided  instantaneous  spectrum  GB{k). 

G7b TFi  =  GA{k)-GB{k)  (2.8) 


The  instantaneous  one-sided  spectrum  for  channel  A  or  B,  GA{k ),  GB(k)  was 
calculated  according  to: 
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Gx(fc)  = 


(2.9) 


SA[k)  for  k  =  0 
2SA{k)  for  1  <  k  <  N/2  -  1 
[  0  for  N/2  <  k  <  N  -  1 

where  N  is  the  number  of  samples  per  ensemble.  For  this  experiment,  N  was 
equal  to  1024.  5x(*)  is  the  two-sided  instantaneous  spectrum  for  channel  A 
which  is  determined  from  the  recorded  time  record,  a(n)  according  to: 

Sa(k)  =  7[u{n)  •  a(n)]  (2.10) 


where  7  is  the  forward  discrete  fourier  transform  and  w(n)  is  the  time  record 
weighting  function.  In  this  experiment,  a  Hanning  window  weighting  function 
was  utilized. 

The  forward  discrete  fourier  transform,  7,  is  defined  as: 

N-l 

=  1/N  a(n)  ■  cxp(~j2nkn/ N)  (2.11) 

n=0 

where  n  equals  the  time  index  times  the  sampling  interval,  nAf  and  k  equals  the 
frequency  index  times  the  frequency  resolution,  fcA/.  The  frequency  resolution 
for  this  experiment  was  2  Hz  for  an  analysis  bandwidth  of  1600  Hz. 

2.5.2  Measurement  techniques 

Two  techniques  can  be  used  to  measure  sound  intensity  using  an  intensity 
probe.  They  are: 

1.  continuous  spatial  averaging 

2.  discrete  point  averaging 

Continuous  spatial  averaging  consists  of  sweeping  the  intensity  probe  at  a 
controlled  rate  over  each  surface  while  continuously  sampling  and  averaging 
data.  This  technique  produces  the  most  accurate  spatially  averaged  results  but 
suffers  slightly  in  repeatability.  Several  repeated  runs,  often  four  or  more,  were 
required  to  obtain  the  required  repeatability  of  total  sound  power  (within  0.5 
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Figure  1:  Louver  Ll  sound  power  spectra  repeatability.  Shown  are  four  succes¬ 
sive  spectral  measurements  of  sound  power. 

dB  on  at  least  three  runs).  Figures  1,  2  and  3  are  the  results  of  three  such 
repeatability  experiments  with  0.5  dB  or  less  variation  in  total  sound  power. 
Shown  is  the  sound  power  spectral  level  from  a  rectangular  orifice,  4  -  1/2"  by 
1  -  2/3",  and  an  analysis  bandwidth  of  2  Hz. 

Discrete  point  averaging  consists  of  dividing  your  control  surface  into  numer¬ 
ous  discrete  areas  and  measuring  the  sound  intensity  at  the  center  of  each  dis¬ 
crete  area  (point)  and  summing  over  all  areas.  This  technique  offers  maximum 
repeatability  but  at  reduced  spatially-averaged  accuracy.  The  two  techniques 
were  both  utilized  initially  and  compared.  The  repeatability  of  the  total  radi¬ 
ated  power  using  the  discrete  point  technique,  with  45  points,  was  0.3  -  0.5  dB 
between  successive  runs.  These  repeatability  values  were  selected  as  the  criteria 
for  the  continuous  spatial  averaging  technique.  The  results  showed  that  the  two 
techniques  gave  similar  sound  power  spectrums  and  similar  total  sound  power 
levels  within  an  accuracy  of  0.3  dB.  For  this  reason,  the  continuous  spatial  aver¬ 
aging  technique  was  selected  due  to  its  increased  accuracy  and  reduced  sampling 
time. 

The  control  surface  was  constructed  of  1/4"  balsa  wood  enclosing  a  total 


9 


ico  403  :cc  eoc 

r>M>ttttr  (i>) 


700  «0O 


Figure  2:  Orifice  01  sound  power  spectra  repeatability.  Shown  are  three  suc¬ 
cessive  spectral  measurements  of  sound  power. 
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Figure  3:  Orifice  01  sound  power  spectra  repeatability.  Shown  are  five  successive 
spectral  measurements  of  sound  power. 
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Figure  4:  Phase  difference  between  two  microphones  and  their  preamplifiers  I 

(channels  A  and  B) 

surface  area  of  0.542  square  meters.  Intensity  spectra  over  the  frequency  range 
of  52-1600  Hz  were  taken  and  recorded  on  a  B  k  K  2032  Spectrum  Analyzer. 

Microphone  calibration  was  performed  using  a  pistonphone  and  daily  tempera¬ 
ture  and  barometer  measurements  were  taken  to  calculate  the  air  density.  To 
ensure  that  the  two  microphone  were  indeed  phase  matched,  their  phase  was 
measured  in  a  random  noise  field  and  compared.  The  microphones  were  placed 
at  the  open  end  of  a  long  cylinder,  equal  distances  from  a  speaker  placed  at  the 
opposite  end  of  the  cylinder.  This  apparatus  was  similar  in  design  to  a  calibra¬ 
tion  pistonphone,  type  4200.  Figure  4  shows  the  phase  difference  measured  over 
a  frequency  range  of  2560  Hz  and  clearly  demonstrates  their  compatibility  over 
the  frequency  band  of  interest,  52  -  1600  Hz. 


I, 
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In  Section  5.2,  the  errors  in  dB  resulting  from  these  phase  mismatches  are 
calculated  and  shown  to  be  on  the  same  order  of  magnitude  as  the  inherent 
statistical  fluctuations. 

2.5.3  Power  spectrum  contamination 

Data  below  52  Hz  were  ignored  due  to  contamination  by  the  1,0,0  room 
mode  (38  Hz)  of  the  semi-anechoic  chamber.  Also  at  this  low  frequenc ,  the 
effectiveness  of  the  anechoic  treatment  was  reduced.  The  peak  at  462  Hz  and  its 
harmonic  at  924  Hz  which  is  evident  in  all  the  spectra  was  due  to  the  resonance 
of  the  cross-duct  mode  of  the  test  section  and  not  at  all  associated  with  the 
orifice  source.  Fortunately,  the  sound  level  and  narrow  bandwidth  associated 
with  this  modal  resonance  contributes  insignificantly  to  the  total  sound  power 
caused  by  the  orifice  source.  Frequencies  above  1600  Hz  contributed  nothing  to 
the  total  sound  power  and  the  spectral  shape  followed  the  same  trend  as  from 
1200  -  1600  Hz. 


3.  THEORY 


3.1  Introduction 


The  develop.nent  of  a  theoretical  model  to  determine  the  scaling  laws  for 
orifice/louver  flow  noise  was  first  accomplished  in  1972  by  Ffowcs  Williams  [8]. 
He  treated  the  model  problem  of  flow  noise  generated  by  a  body  of  turbulence 
near  an  infinite,  rigid,  thin  screen  with  circular  orifices.  His  assumptions  were: 
a)  the  size  of  the  orifice  is  much  smaller  than  the  distance  between  them,  b) 
variations  in  pressure  within  individual  apertures  are  neglected.  While  these 
assumptions  are  rather  restrictive  in  their  application  to  many  industrial  louvers, 
the  theory  does  provide  a  fundamental  result:  that  the  sound  radiated  by  orifices 
and  louvers  can  be  directly  related  to  the  mass  flow  driven  through  the  aperture 
by  the  aerodynamic  field. 

This  fundamental  result  was  used  by  Nelson  [9]  to  develop  a  model  for  flow 
noise  generated  over  practical  perforated  screens  with  circular  orifices.  Due  to 
the  normal  component  of  the  turbulent  fluctuations  above  the  louver,  flow  is 
driven  through  the  orifice  resulting  in  a  fluctuating  volume  flow  through  the 
upper  surface.  This  fluctuating  volume  flow  constitutes  an  equivalent  monopole 
source  at  the  upper  surface.  Since  the  volume  flow  into  the  upper  surface  of 
the  orifice  creates  an  equivalent  outflow  through  the  lower  surface,  there  is  an 
equal  and  opposite  equivalent  monopole  at  the  lower  surface.  The  net  field 
which  radiates  out  of  the  louver  and  into  the  closed  acoustic  space  is  thus  the 
combination  of  the  lower  monopole  source  and  the  contribution  from  the  upper 
monopole  source  which  is  dependent  on  the  degree  of  acoustical  transparency  of 
the  louver.  If  the  louver  is  acoustically  opaque,  then  only  the  lower  monopole 
source  will  radiate  into  the  acoustic  space.  If  the  louver  is  acoustically  trans¬ 
parent,  then  the  upper  monopole  source  will  radiate  into  the  space,  having  a 
cancelling  effect  on  the  field  of  the  lower  monopole  source.  This  would  result  in 
a  net  radiated  field  which  exhibits  properties  of  a  dipole  source. 

Nelson’s  equivalent  source  model  gave  an  excellent  description  of  the  flow 
noise  produced  by  perforated  facing  materials  used  in  practice  even  though  he 
assumed  a  rather  simplistic  form  for  the  driving  aerodynamic  velocity.  His 
assumptions  were:  a)  the  wavelength  for  the  sound  is  much  larger  than  the  orifice 
radius,  b)  there  are  many  apertures  per  acoustic  wavelength  and  no  interactions 
of  the  flow  occurred  between  adjacent  apertures  which  can  modify  the  aperture’s 
conductivity.  In  the  sections  that  follow,  the  theory  developed  by  Nelson  and 
Ffowcs  Williams  will  be  utilized  to  model  the  flow  noise  created  by  turbulent  flow 


13 


Figure  5:  Coordinate  system  for  analysis 

over  orifices  and  louvers.  Once  the  equivalent  source  model  has  been  developed, 
it  will  be  used  to  collapse  the  experimental  sound  power  spectra  thus  providing 
a  verification  to  the  theory’s  validity.  The  coordinated  system  used  for  this 
analysis  is  shown  in  Figure  5. 

3.2  The  Equivalent  Source  Model  for  Orifices  and  Louvers 

Ffowcs  Williams  [8]  and  Leppington  and  Levine  [10]  developed  an  expression 
for  the  reflection  coefficient,  R,  for  a  plane  wave  of  single  frequency,  w,  incident 
at  a  polar  angle,  0,  on  an  infinite  louver.  For  low  frequencies  and  orifice  spacings 
that  are  large  compared  to  the  orifice  size,  the  reflection  coefficient,  R,  is 

R  =  1/(1  +  4  aiN/kcosO)  (3.1) 


where  a  is  the  orifice  radius,  N  the  number  of  orifices  per  unit  area  of  louver 
and  k  is  the  acoustic  wavenumber,  k  =  ufCo. 

The  field  produced  by  the  lower  monopole  at  y*  can  be  determined  using  the 
principle  of  reciprocity.  The  field  at  y*  due  to  a  distant  source  at  x*,  is  given  by 
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Figure  6:  The  infinite  screen  reflection  coefficient  as  a  function  of  the  ratio 
(AaN/kcosO).  From  [9,  Nelson,  P.A.,  1982]. 

p,(  1  +  R)  where  p,  is  the  incident  plane  wave  striking  the  lower  surface  of  the 
louver.  Thus  one  can  write 

P{Vk)  =  (1  +  R)Q(xk)exp(ik  1  yk  -  xk  \)/An  \  yk  -  xk  \  (3.2) 


where  Q  is  the  strength  of  the  distance  monopole  at  xk.  By  reciprocity,  the  field 
at  xk  due  to  the  monopole  at  yk  is  given  by  interchanging  the  source  and  receiver 
position  vector  in  equation  3.2.  Equation  3.2  can  be  generalized  to  relate  the 
spectral  density  of  the  acoustic  pressure  field,  Sp(xk,u)  to  the  spectral  density 
of  the  source  Sq( y*,w)  such  that 

S„(z*,W)  =|  1  +  R  |*  -SQ{Ykyuj)l  16*V  (3.3) 


where  r 1  =|  xk-yk  |*  .  Figure  6  is  a  plot  of  |  R  |l  against  the  ratio  (AaN /k  cos  9). 
From  this  plot,  it  can  be  seen  that, 

j  J*=  1  for  (4  aN/k  cos  9)  <  0.3 
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This  occurs  for  open  area  ratios  on  the  order  of  0.01  (see  reference  11).  For  all 
the  orifice  and  louvers  tested  in  this  study,  the  open  area  ratio  is  less  than  0.005 
so  that  j  R  |*=  1.  To  verify  this  result  for  the  reflection  coefficient,  since  the 
assumption  of  many  apertures  per  acoustic  wavelength  is  not  valid  for  the  single 
orifice,  a  separate  analysis  was  performed.  Here,  in  two  dimensions,  the  plane 
wave  of  single  frequency, u,  is  incident  normal  on  an  infinite,  rigid  surface  with 
a  single  opening,  of  dimension,  2b.  Again  the  principle  of  reciprocity  is  used  to 
predict  the  radiated  pressure  from  the  single  orifice.  Using  matched  asymptotic 
expansions,  for  6/ A  <<  1,  the  lowest  order  solution  for  the  reflection  coefficient 
is  again,  |  R  |J=  1.  The  details  of  this  analysis  are  shown  in  Appendix  1. 

Thus, 


Sp(zt,u)  =  S<3(yk,w)/47rV 


(3.4) 


The  instantaneous  source  strength  <j(yk,t)  can  be  written  in  terms  of  the 
weighted  integral  of  the  unsteady  normal  component  of  velocity  over  the  area  of 
the  aperture.  To  take  advantage  of  the  homogeneous  and  stationary  properties 
of  the  velocity  in  the  plane  of  the  wall,  we  can  write 

?(y*.0  =  -*wp  f  f  v{yk,t)K(sk~yk)dA(sk)  (3.5) 

J  OO  j  OO 


where  yk  is  the  vector  position  of  some  reference  point  on  the  aperture  surface 
and  K(sk  -  yk )  is  the  weighting  function,  i.e.,  the  function  which  converts  the 
integral  over  the  aperture  area,  sk  =  ,si,sj,  to  an  integral  over  all  space,  defined 
as 


K(sk  —  yk)  =  1  for  sk  -  yk  inside  of  the  aperture  area,  A 

=  0  for  sk  -  yk  outside  A 

Following  Corcos  [12]  the  spectral  density  of  the  source  strength,  S<}(yk,w) 
can  be  defined  by 

Sq(y*.w)  =  p*w2  f  [  Svv(zk  -  sk,u)K[zk)K{sk)dA{zk)dA{sk)  (3.6) 

J  OO  j  OO 
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where  zk  is  a  dummy  space  vector.  Define  ek  by 


« k  =  Zk~  ■** 


thus 

S<){yk,u)  =  pV  [  [  Svv[ek,u)K{sk)K{sk  +  tk)dA(sk)dA(ek) 

J OO  y  oo 

=  pJwJ  /  5uu(£k,w)^(£t)di4(fi)  (3.66) 

J  OO 

where  the  function,  #(£*),  defined  as  the  orifice  response  function,  is  determined 
wholly  by  the  geometry  of  the  aperture  and  can  be  evaluated  as 

/oo 

K(sk  +  £4)di4(s*)  =  1  x  area  of  overlap 

=  (L-  |  f3  |)(6~  |  |)  for  |  £3  |<  L,\  £j  |<  6  (3.6c) 

=  0  elsewhere 

for  the  rectangular  shaped  orifice  shown  below 


Finally,  we  assume  the  cross-spectral  density,  S„u(£*,u;J,  to  have  the  product 
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form 


w)  =  Sv{u)A{Kt1)B{Kis)eiK^.  (3.7) 


where  S„(u/)  is  the  spectral  density  of  the  normal  velocity  fluctuations,  K  is 
the  convective  wavenumber  given  by  K  =  u >jUe  where  Ue  is  the  convection 
speed  of  the  boundary  layer  eddies  and  A[K(i )  and  B(Kt3)  are  the  amplitude 
functions  for  the  normalized  longitudinal  and  lateral  cross-spectra,  respectively. 
This  form  is  similar  to  the  expression  given  by  Corcos  [  12]  for  the  cross  spectral 
density  of  the  boundary  layer  pressure  fluctuations.  The  longitudinal  cross¬ 
spectra  of  the  eddies  as  they  convect  in  the  yi  direction  also  contains  the  phase 
term  txp{iKt\).  A[Ke\ )  is  plotted  in  Figure  7  as  a  function  of  its  argument 
for  various  ratios  of  separation  distance  to  momentum  thickness  ( ry/6 *)  for  wall 
pressure.  Since  the  wall  pressure  is  related  to  the  volume  integral  of  the  normal 
component  of  the  velocity,  wall  pressure  cross-spectra  data  can  be  used  to  predict 
the  longitudinal  and  the  lateral  cross-spectra  of  the  normal  component  of  particle 
velocity.  Comparing  the  wall  pressure  correlation  length  scales  from  1 1 2]  to  the 
velocity  correlation  length  scales  from  [l7j,  we  find  that  both  the  longitudin'al 
and  lateral  length  scales  agree  to  within  an  order  of  magnitude,  thus  supporting 
the  previous  statement. 

For  the  various  orifice  and  louver  samples  tested,  at  a  maximum  value  for 
ri/6*  equal  to  1,  A(Kex)  varies  from  a  maximum  value  of  1  to  a  minimum 
value  of  0.6  for  orifice  02  at  U  =  20 m/s  and  /  =  250  Hz  (the  freque  y  at 
which  the  spectral  level  is  lOdB  down  from  the  peak  level).  Throughout  most 
of  the  samples  tested,  A(Kt\)  =  1  over  the  frequency  band  of  interest.  Thus 
the  longitudinal  correlation  lengths  were  typically  several  times  the  length  of 
the  orifice  aperture.  Thus  the  assumption  is  made  that  no  decay  in  the  eddy 
structure  occurs  in  the  longitudinal  (yt)  direction. 

The  normalized  lateral  cross-spectra  in  the  y3  direction  is  given  by  the  term 
B(Ke3)  which  is  assumed  unity  for  no  separation  and  to  decay  very  rapidly  with 
increased  separation.  To  justify  this  assumption,  B[Kt 3)  is  plotted  in  Figure 
8  as  a  function  of  its  argument  for  various  ratios  of  r3/6*.  From  this  figure, 
the  typical  correlation  length,  l,  is  calculated  to  be  on  the  order  of  2  cm.  The 
ratio  of  the  lateral  correlation  length,  l,  to  the  aperture  length,  L,  is  0.18  thus 
supporting  the  initial  assumption  of  rapid  lateral  decay.  The  correlation  lengths 
are  determined  at  the  point  where  the  amplitude  function  equals  0.4  and  at  the 
frequency  corresponding  to  the  peak  power  spectral  level. 
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Figure  7:  Amplitude  function  for  longitudinal  cross-spectral  density  of  wall  pres¬ 
sure.  From  [16,  BBN  report  No.  4110366,  1966]  I 
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Figure  8:  Amplitude  function  for  lateral  cross-spectral  density  of  wall  pressure. 
From  [16,  BBN  Report  No.  4110366,  1966] 


Substituting  equation  3.7  into  equation  3.6b  and  evaluating  the  integral  over 
the  rectangular  shaped  orifice,  assuming  Sv(yktu)  to  be  uniform  over  the  orifice 
area  gives: 


S<j(y*,u;)  =  p2w*S„(y*,w)  •  Li  •  (2//tT2)(l  -  cos{Kb ))  (3.8) 


where  L  is  the  orifice  length  in  the  ys  direction.  See  Appendix  (3)  for  the 
details  of  the  integration  for  the  rectangular  shaped  orifices  and  reference  [9]  for 
the  results  for  circular  shaped  orifices.  Equation  3.8  states  that  the  monopole 
spectral  density  is  directly  related  to  the  normal  velocity  spectral  density. 

3.3  Scaling  Laws  for  the  Radiated  Sound  Power 

Having  determined  the  form  for  the  source  strength,  we  can  now  develop 
the  scaling  laws  for  the  radiated  sound  power  from  various  rectangular  shaped 
configurations.  Assuming  that  each  orifice  radiates  sound  incoherently,  the  to¬ 
tal  source  strength  for  N  apertures  is  simply  N  times  the  individual  source 
strength  for  each  aperture  defined  by  equation  3.8.  Substituting  equation  3.8 
into  equation  3.4  gives: 

Sp( i*,w)  =  Np2Ll  •  Su(yt,w)(l  -  cos{Kb))U2 /2jtV  (3.9) 


The  assumption  is  made  that  the  convection  velocity,  Uc,  is  equal  to  one  half 
the  freestream  velocity,  Ua 0.  Define  the  turbulent  intensity,  /,  and  the  Strouhal 
number,  St(£)  and  St(b)  by: 

(/’)a«  =  1/Ul  ■  j  Sv{yk, u)du>  (3.10) 


St(t)  =  ft/U00,St(b)=  fb/U,,  (3.11) 


Finally,  integrating  equation  3.9  over  the  bandwidth  Aw  with  the  above 
expressions  substituted,  yields 
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(3.12) 


p’(**.')U.  =  Np'L  ■  (1  -  coi(«-6))SI(<)(/s)i„t/i/2/^rI 


For  spherical  waves,  the  sound  power,  P,  can  be  determined  by: 

C PWLu  =  NpL{  1  -  cosiKbfiStmPUMfCo  (3.13) 


Equation  3.13  is  the  final  result  for  the  flow  noise  created  by  turbulent  flow  over 
orifices  and  louvers  with  rectangular  shape  as  described  in  Section  (4.2).  The 
controlling  parameters  for  scaling  sound  power  are  thus: 

1.  Strouhal  number  based  on  the  lateral  correlation  length,  l  and  aperture 
transverse  dimension,  6 

2.  the  normal  component  of  the  Turbulent  Intensity 

3.  the  fifth  power  of  the  freestream  velocity 

4.  the  shape  factor  (1  -  cos(/f6)) 

In  Section  (5.10),  this  model  is  used  to  collapse  the  experimental  sound 
power  data  for  orifices  02  and  03  and  louver  L 1.  It  is  clear  from  the  collapse 
that  equation  3.13  does  indeed  identify  the  predominant  controlling  parameters 
for  rectangular  shaped  orifice  and  louver  flow  noise.  Note,  for  circular  shaped 
apertures,  the  above  scaling  factors  remain  the  same  except  the  shape  factor  is 
replaced  by  the  function,  H\{Kd)  where  Hl  is  a  Struve  function  of  the  first  kind 
and  d  is  the  aperture  diameter.  See  reference  9  for  the  details  in  calculating 
equation  3.8  for  a  circular  aperture. 
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4.  EXPERIMENTAL  APPARATUS 


4.1  Semi-anechoic  Wind  Tunnel 

The  experiments  were  carried  out  using  the  low  turbulence,  sub  sonic  wind 
tunnel  of  the  Acoustics  and  Vibration  Laboratory  at  MIT.  The  wind  tunnel  de¬ 
tails  are  shown  in  Figure  9  with  the  basic  construction  characteristics  described 
by  Hanson  [13]  and  modifications  described  by  Shapiro  [14].  The  wind  tunnel  is 
of  open  circuit  design  with  the  test  section  enclosed  in  a  blockhouse.  The  tunnel 
consists  of  a  set  of  flow  straighteners,  screens,  a  settling  chamber  and  a  20:1  con¬ 
traction  leading  to  a  square  test  section  38cm  by  38cm.  The  blockhouse  walls, 
floor  and  ceiling  were  treated  with  urethane  foam  to  provide  a  semi-anechoic 
chamber  within  which  acoustic  intensity  measurements  were  taken.  Figure  10 
details  a  cut-away  front  and  top  view  of  the  chamber  and  test  section.  The 
acoustic  treatment  consisted  of  a  base  of  10cm  foam,  then  a  set  of  10cm  foam 
cubes  spread  randomly  and  finally  a  set  of  2cm  jam  cubes  also  spread  randomly. 
The  effectiveness  of  the  coating  was  determined  by  [14]  using  a  horn  driver  and 
measuring  the  sound  field  decay  with  distance.  The  anechoic  treatment  was 
effective  between  200  Hz  and  2000  Hz  and  only  partially  effective  at  100  Hz.  ' 

4.2  Orifice/Louver  Samples 

Nine  separate  orifice/louver  samples  were  tested.  All  samples  were  made  of 
plexiglass  and  screwed  into  the  sample  mounting  window  which  was  mounted 
flush  into  the  side  to  the  test  section  as  shown  in  Figure  10.  The  separate 
samples  tested  were: 

1.  O-l:  a  4  -  1/2"  (height)  by  1  -  2/3"  (width)  rectangular  orifice,  3/4"  thick 
with  straight  sided  walls 

2.  0-2:  the  above  made  of  1/4"  thick  plexiglass 

3.  0-3:  a  4  —  7/16M  by  27/64"  rectangular  orifice,  1/4"  thick  with  straight 
sided  walls 

4.  0-2A:  0-2  with  40  degree  slanted  walls,  i.e,  both  walls  form  an  angle  of 
40  degrees  towards  upstream  from  the  perpendicular  to  the  flow. 
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Figure  9:  Wind  Tunnel  Facility 

6.  L-l:  4  rectangular  orifices,  each  4-7/16"  by  27/64",  separated  by  1/4",  1/4" 
thick 

7.  L-2:  54  circular  orifices,  each  27/64"  in  diameter,  separated  by  0.6328" 
between  centers,  1/4"  thick 

8.  L-1A:  L-l  with  40  degree  slanted  walls 

9.  L-2A:  L-2  with  40  degree  slanted  walls. 

The  louver  configurations  were  designed  to  maintain  the  same  total  open 
area  as  the  single  orifices,  0-1  and  0-2.  The  40  degree  slanted  walls  provided  a 
blunt  orifice  trailing  edge  and  subsequently  a  sharp  orifice  leading  edge,  as  seen 
from  the  flow. 

4.3  Test  Section 

The  test  section  was  made  from  3/4"  finished  plywood,  38  cm  by  38  cm 
square  in  cross  section  and  was  located  on  the  centerline  of  the  chamber.  The 
test  section  was  sealed  and  the  ends  were  further  covered  by  2"  of  urethane  foam. 
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Figure  10:  Experimental  set-up  within  the  wind  tunnel,  front  and  top  views 
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The  background  power  levels,  measured  with  the  orifice  covered,  were  sufficiently 
low  compared  to  the  power  levels  with  the  orifice  exposed  (see  Section  5.3)  to 
ensure  signal  to  noise  ratios  of  at  least  10  dB  over  the  frequency  range  of  interest 
and  at  all  operating  speeds. 

4.4  Turbulent  Intensity  Measurements 

Turbulent  Intensity  is  defined  by  the  relationship: 

77  =  20log(v^/t/0o  (4.1) 


Turbulent  fluctuations  were  measured  using  a  Dantec  Constant  Temperature 
Anemometer  with  a  5  micron  diameter  tungsten  wire.  The  output  of  the 
anemometer,  being  the  total  voltage  E,  was  inputted  directly  to  the  B  &  K  2032 
Spectrum  Analyzer.  Here  the  AC-coupled  autospectrum,  Ge(uj),  was  measured. 
The  mean-squared  fluctuating  velocity  autospectrum,  Gu(w)  was  determined 
from: 


Gu[u)  =  (\H  r)G*M 


(4.2) 


where  H  is  the  transfer  function.  See  Appendix  2  for  the  details  in  determining 
H.  Turbulent  Intensity  spectra  were  measured  at  all  four  flow  speeds  for  orifice 
0-2  and  louvers  L  -  1  and  L  —  1  A.  The  position  of  the  U  probe  was: 

1.  two-tenths  of  an  inch  into  the  boundary  layer  directly  above  the  orifice  at 
both  the  upstream  and  downstream  edge.  For  the  louvers,  the  streamwise 
measurements  were  taken  in  the  center  of  the  furthermost  upstream  and 
downstream  orifice, 

and 

2.  half-way  into  the  throat  of  the  orifice  and  at  the  same  streamwise  locations 
as  described  above. 
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5.  EXPERIMENTAL  RESULTS 


5.1  Introduction 


Figure  11  shows  the  geometry  of  the  problem.  Air  flows  at  free  speed  Uoo  over 
an  orifice  or  louver  which  has  been  inserted  flush  into  a  long  flat  walled  duct.  The 
orifice  was  either  rectangular,  with  streamwise  dimension  b  and  length  L  or  cir¬ 
cular  with  diameter  d.  The  acoustic  radiated  power  from  the  orifice/louver  was 
measured  at  free  stream  velocities  of  20,  30,  40  and  50  meters  per  second.  The 
acoustic  power  was  determined  by  sweeping  the  intensity  probe  over  a  control 
surface  (a  five  sided  boxed  frame  surrounding  the  orifice)  then  multiplying  the 
resulting  spatially  and  temporally  averaged  intensity  by  the  total  surface  area 
of  the  control  surface.  Five  orifices  and  four  louvers  were  tested  as  described  in 
Section  4.2. 

To  determine  the  contribution  made  solely  by  the  orifice/louver  source,  two 
additional  experiments  were  conducted.  One  experiment  consisted  of  measuring 
the  acoustic  power  radiated  by  the  duct  when  the  orifice  was  covered.  The 
second  experiment  consisted  of  measuring  the  power  flow  through  the  orifice,  in 
the  absence  of  free  stream  flow,  due  to  an  artificially  created  acoustic  field  in 
the  duct.  This  field  was  created  using  a  speaker  in  the  duct  and  was  equal  in 
magnitude  to  that  caused  by  the  turbulent  boundary  layer  at  U  =  50  m/s.  This 
experiment  determined  the  “open-window”  contribution  to  the  overall  sound 
power.  Finally,  directivity  measurements  were  taken  for  orifice  0-2  and  louver 
L  -  2 A  at  U  =  25  m/s  and  the  respective  polar  plots  produced. 

5.2  Reactivity  of  the  Sound  Field  and  Sound  Intensity  Error  Limits 

The  reactivity  of  the  sound  field  in  the  blockhouse  (the  closed  acoustic  space) 
was  measured  for  various  orifice  and  louver  samples  (see  Figures  12  and  13).  The 
Table  below  lists  the  measured  difference  between  the  total  sound  pressure  level, 
Lp  and  the  total  sound  intensity  level,  £,<,  for  the  orifice  02A  and  the  louver 
L2A  at  a  radial  distance  of  2  feet. 


Sample 

Reactivity  ( Lp  -  Li) 

02A 

0.8 

MEM 

1.0 

BUI 

L2A 

0.9 

msM 

1.0 

Table  1:  Measured  Sound  Field  Reactivity 
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Figure  11:  Schematic  of  orifice  and  probe  geometry.  Air  flows  at  freestream 
velocity,  U,  over  an  aperture  of  streamwise  dimension,  b.  Sound  intensity  is 
measured  by  probe  (microphone  spacing,  Ar)  while  continuously  sweeping  con¬ 
trol  surface. 


28 


tU - - - - - - - - - - - - - - - 1 

□  D.  2H  O.  ih  0.  6h  □.  ah  1  .  Oh  )  .  2h  ].4h  1  .  6h 

Frequency  (Hi) 


Figure  12:  Sound  field  reactivity  for  orifice  02A.  Shown  are  the  sound  intensity 
and  sound  pressure  spectral  levels.  The  difference  between  the  two  is  defined  as 
the  sound  field  reactivity. 
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Figure  13:  Louver  L2A  sound  intensity  spectral  levels  versus  sound  pressure 
spectral  levels  ( Lfvs.Lp ).  The  difference  between  Li  and  Lp  is  defined  as  the 
sound  field  reactivity. 
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From  Figure  4,  the  measured  phase-mismatch  for  the  sound  intensity  probe 


Frequency  (Hz) 

4>m{degrees) 

63 

0.8 

500 

0.2 

1000 

0.5 

2000 

0.9 

The  maximum  phase-mismatch  between  channels  A  and  B  for  the  Analyzer 
is  0.3  degrees.  With  a  spacer  value  of  50mm  and  the  above  sound  field  reactivity, 
the  phase  reactivity  nomogram  of  reference  6  gives  the  following  value  for  the 
maximum  phase  of  the  sound  field,  between  the  two  microphone  positions. 


Frequency  (Hz) 

<j>m(degrees) 

63 

2 

500 

16 

1000 

25 

2000 

60 

Therefore  the  limits  of  the  phase-mismatch  are: 


Frequency(Hz) 

4>f  ~  <t> a  -  4>m 
(degrees) 

4>J  +  <t>a  +  4>m. 
(degrees) 

63 

0.9 

3.1 

500 

15.5 

16.5 

1000 

24.2 

25.8 

2000 

58.8 

61.2 

The  above  total  phase-mismatch  values  are  the  limits  of  error  in  the  measured 
phase  at  the  frequency  listed.  The  corresponding  limits  in  dB  for  the  measured 
intensity  from  the  nomogram  are: 


Frequency  (Hz) 

Error  limits  (dB) 

63 

+  2.5 

500 

+  0.5 

1000 

+  0.25 

2000 

+  0.125 

J 


31 


The  corresponding  statistical  standard  deviation,  f ,  of  the  measured  spectral 
levels  about  the  mean  value  is: 


€  =  4.34  /jBTa,{dB) 


(5.1) 


where  B  is  the  filter  bandwidth  and  Ta  the  total  time  record  length.  Substituting 
the  appropriate  values  for  B  (2  Hz)  and  Ta  (150  sec)  gives  e  =  0.25  dB.  Thus 
the  limits  of  error  in  measuring  the  sound  intensity  are  approximately  equal  to 
the  inherent  statistical  fluctuations  (for  frequencies  above  100  Hz). 

5.3  Length  and  Frequency  Scaling 

The  three  significant  length  scales  for  this  experiment  were: 

1.  the  streamwise  orifice  dimension,  6 

2.  the  turbulent  boundary  layer  thickness,  <5,  above  the  orifice,  and 

3.  the  lateral  correlation  length,  t. 

Table  2  provides  the  applicable  values  for  6,  6  and  l  used  in  this  exper¬ 
iment.  The  lateral  correlation  length  i  was  inferred  from  existing  transverse 
cross-spectra  of  wall  pressure  beneath  turbulent  boundary  layers  (see  Appendix 
3). 


Orifice/ 

Louver 

b(cm) 

U(m/s) 

l  (cm) 

6  (cm) 

6/b 

for  0-2 

6/b 

for  0-3 

0-1 

4.23 

20 

2.4 

4.19 

0.9» 

3.95 

0-2./0-2A 

4.23 

30 

1.8 

3.86 

0.91 

3.64 

O-3./0-3A 

1.07 

40 

0.9 

3.65 

0.86 

3.44 

L-1,/L-1A 

1.07 

50 

1.4 

3.49 

0.83 

3.29 

L-2./L-2A 

1.07 

50 

2.1 

3.49 

Full  Scale 

2.50 

16.5 

4.70 

1.88 

Table  2:  Experimental  Length  Scales 

Note:  Full  scale  assumes  a  louver  positioned  15  ft  from  the  leading  edge  of  a  flat  surface 
in  seawater  at  58  F.  The  louver  consists  of  multiple  rectangular  shaped  openings,  each 
of  transverse  dimension,  6  =  2.5  cm.  From  Figure  16  at  6/b  =  1.88,  the  Strouhal  #  = 
0.19.  Therefore  the  full  scale  frequency,  /e,  is: 


fe  -  (0.19) (16.5)/(0.025)  =  125  Hz 
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The  significant  frequency  scales  were  the  non-dimensional  strouhal  numbers 
St(b )  and  St(£)  defined  by: 


S*(6)  =  feb/U » 

St{l)  =  ul/Uc  =  Anfcl/Ugo  (5.2) 


where  fe  is  the  center  frequency  of  the  1/3  octave  band  with  the  peak  sound 
power  level.  Figure  14  is  a  graph  of  fe  versus  for  orifice  0  —  1  at  a  constant 
b.  The  data  points  reflect  several  separate  experiments  measuring  the  sound 
pressure  level  at  radial  distances  of  9,  18,  24  and  36  inches  from  the  orifice. 
Figure  15  plots  fc  versus  {/*,  for  orifices  0-1, 0-2, 0-3  and  louver  L  —  1  and 
L  -  2.  Clearly  there  exists  a  linear  relationship  (or  nearly  linear)  between  fc 
and  Ugo  for  a  particular  transverse  dimension  b. 

Table  3  lists  the  non-dimensional  Strouhal  numbers  for  the  above  ori¬ 
fices/louver  as  well  as  the  non-dimensional  length  scale,  6/b  for  U  —  20  m/s. 
The  strouhal  number,  St  =  fb/U,  has  been  defined  to  equal  that  recommended 
by  DeMetz  and  Farabee  [5]  to  compare  results.  Figure  16  is  the  graph  of  the 
data  in  Table  3  along  with  the  results  of  DeMetz  and  Farabee  [5]  for  TBL  flow 
over  orifices  enclosed  by  solid  cavities. 


Oiifice/ 

Louver 

5t(6)  =  fb/U 
Strouhal  # 

Non-dimensional 
length  scale  6/b 

O-l 

0.17 

1.22 

0-2 

0.23 

1.22 

0-3 

0.12 

4.79 

0.076 

4.79 

0.048 

4.79 

Table  3:  Non-dimensional  Frequency  and  Length  Scales  at  U  =  20  m/s 

The  decreasing  trend  in  the  Strouhal  number  with  increasing  6/b  for  the 
two  cases  agrees  reasonably  well.  Having  a  knowledge  of  the  boundary  layer 
thickness,  one  can  enter  Figure  16  at  the  appropriate  value  for  6/b  and  determine 
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Figure  14:  Frequency /velocity  relationship  for  orifice  01.  Shown  are  four  sepa¬ 
rate  experiments  for  various  radial  distances,  r. 
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Figure  15:  Frequency/velocity  relationships  for  orifices  01,02,03  and  louvers 

Ll,L2 
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Figure  16:  Strouhal  number  versus  boundary  layer/aperture  width  ratio 
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the  Strouhal  number  for  orifice  or  louver  flow  noise.  From  the  Strouhal  number 
and  the  geometry  of  the  aperature(s),  one  can  predict  fe. 

5.4  Structural  Radiation  Due  to  TBL  Excitation 


To  determine  the  contribution  from  structural  radiation  caused  by  the  exci¬ 
tation  from  the  TBL,  an  experiment  was  conducted  at  all  free  stream  velocities 
with  the  orifice/louver  covered.  Identical  sound  power  spectrums  were  recorded 
as  described  in  Section  4.1.  The  resulting  background  power  spectrums  levels, 
Figures  17,  18,  19  and  20,  were  10  to  46  dB  lower  than  with  the  orifice  exposed. 
From  these  figures,  the  major  contribution  from  the  orifice/louver  sources  was 
in  the  range  of  65  to  1000  Hz.  Background  noise  was  within  15  dB  of  the 
orifice/louver  noise  below  50  Hz  and  above  1000  Hz. 

5.5  “Open  Window”  Contribution 


Previous  work  by  Spinka  [15]  determined  the  acoustic  pressure  field  within 
the  test  section  at  various  free  stream  velocities.  To  determine  the  contribution 
from  the  open  window  affect,  a  similar  acoustic  pressure  field  was  generated 
within  the  duct  using  a  speaker  driven  a  B  &  K  random  noise  generator,  type 
1402,  and  a  B  &  K  band  pass  filter  set,  type  1612.  The  resulting  power  which 
radiated  through  the  orifice  was  measured  as  described  in  Section  2.5  and  com¬ 
pared  to  the  power  radiated  by  the  orifice  source.  Table  4  lists  the  resulting 
power  levels  for  both  cases  at  U  =  50  m/s.  The  conclusion  was  that  the  open 
window  affect  results  in  levels  which  were  typically  23  dB  lower  than  the  levels 
due  to  the  orifice  source. 


1/3  OB  center 
frequency 

Open  window 
power  in  1/3  OB  (2) 

Orifice  O-l 
power  in  1/3  OB  (3) 

(3-2) 

400 

53.8 

77.6 

23.8 

500 

48.2 

72.7 

24.5 

630 

45.9 

69.2 

23.3 

800 

43.7 

67.9 

24.2 

1000 

46.2 

66.4 

20.2 

1250 

41.0 

64.5 

23.5 

Table  4:  “Open  window*  power  level*  versus  orifice  O  -1  power  levels  at  U  =  50  m/s 


5.6  Orifice  Source  Power  Spectra  and  Total  Power 

The  measured  sound  power  spectra  versus  the  free  stream  velocity,  U  = 
20,30,40  and  50  m/s,  for  O  —  1,0  —  2, 0  —  3,  L  -  1  and  L  -  2  are  shown  in 

Figures  21,  22,  23,  24,  25.  The  total  sound  power  radiated  by  each  orifice/louver 
source  at  the  above  speeds  is  given  in  Table  5. 
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Figure  21:  Orifice  Ol  sound  power  spectral  level  versus  free  stream  velocity 
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Figure  22:  Orifice  02  sound  power  spectral  levels  versus  freestream  velocity 
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Figure  23:  Orifice  03  sound  power  spectral  levels  versus  freestream  velocity 
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Figure  24:  Louver  Ll  sound  power  spectra  levels  versus  freestream  velocity 
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Figure  25:  Louver  L2  sound  power  spectral  levels  versus  freestream  velocity 


2000. 
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Table  5:  Orifice/louver  total  sound  power  levels  versus  freestream  velocity,t/oo 


Several  observations  are  made  from  reviewing  the  above  data. 

.1  The  effect  of  louvering  on  sound  power 

•  louver,  L-l,  radiated  the  most  acoustic  power,  averaging  10  dB  more  at  all 
speeds  than  the  single  orifice  though  the  total  open  areas  were  the  same 
(see  Figures  26  and  27). 

•  louver,  LI,  showed  strong  resonant  peaks  and  had  the  narrowest  half-power 
bandwidth  of  all  the  samples  tested. 


Figure  26:  The  effect  of  louvering  on  sound  power,  L 1  versus  02  at  U  =  30  m/s 

•  louver,  Ll,  radiated  significantly  more  power  than  louver  L2  in  the  lower 
frequency  bands  but  approximately  equal  power  in  the  middle  to  high 
frequency  bands  (see  Figures  28  and  29) 

•  louver,  L2,  had  the  broadest  spectra  of  all  the  samples  tested  with  equal 
open  area. 

•  louver,  L2’s  total  radiated  power  was  within  2  dB  of  orifice  02’s  power  in 
all  speed  ranges  but  the  spectral  shape  was  significantly  broader.  As  shown 
in  Figures  30  and  31,  louver  L2  radiated  significantly  more  power  than  02 
in  the  middle  to  higher  frequency  bands  (10  dB  more  on  average)  but 
slightly  less  (approximately  1-2  dB)  power  in  the  lower  frequency  bands. 

Thus  the  effect  of  louvering,  simply  by  dividing  the  total  open  area  into  more 
numerous,  smaller  areas  does  not  by  itself  reduce  the  total  radiated  flow  noise 
and  will  most  probably  result  in  increased  radiated  sound  power  in  the  higher 
frequency  bands. 


2008. 
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Figure  28:  Louver  geometry  effect  on  sound  power,  L 1  versus  L2  at  U  —  30  m/s 
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Figure  29:  Louver  geometry  effect  on  sound  power,  L 1  versus  L2  at  U  —  50  m/ s 
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Figure  32:  Orifice  depth  effect  on  sound  power,  Ol  versus  02  at  U  —  30  and 
50  m/s 

5.6.2  Orifice  throat  thickness  effect  on  radiated  power 

Figure  32  shows  the  orifice  throat  depth  effect  on  sound  power  by  comparing 
orifice  Ol  to  orifice  02.  Orifice  02  radiated  more  power  (approximately  3dB) 
in  the  higher  frequency  bands  but  maintained  approximately  1-2  dB  less  power 
in  the  more  dominant  lower  frequency  bands. 

Thus  the  throat  thickness  was  more  predominant  in  effecting  the  spectral 
shape  than  reducing  the  total  radiated  power.  The  effect  of  reducing  the  orifice 
throat  thickness  appears  to  broaden  the  spectral  shape  while  maintaining  the 
same  total  radiated  sound  power. 

5.6.3  Streamwise  dimension,  b,  effect  on  sound  power 

The  spectral  shape  of  orifice  03  had  the  widest  half-power  bandwidth.  The 
effect  of  reducing  the  streamwise  dimension,  b,  was  to  broaden  considerably  the 
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power  spectra  while  reducing  only  moderately  the  total  sound  power.  From 
section  3,  the  controlling  shape  factor  was  (l-cos(Kb)).  Thus  the  predicted 
result  from  reducing  b  by  four  fold  would  be  a  subsequent  increase  by  four  fold 
in  the  power  spectra  bandwidth.  Reviewing  again  the  power  spectra  for  orifices 
02  and  03  (Figures  22  and  23),  we  see  that  indeed  the  power  bandwidth  for  03 
is  approximately  four  times  that  of  02  (defining  the  bandwidth  to  be  between 
the  one-quarter  peak  spectral  power  points). 

5.7  Speed/Power  Laws 

The  speed/power  law  is  of  the  form: 


P  =  CUn 


where  C  is  a  value  which  is  independent  of  speed.  The  data  in  Table  5  was  uti¬ 
lized  to  determine  the  total  radiated  power  dependence  on  speed.  The  equations 
used  to  determine  the  value  of  the  speed  exponent,  n,  were: 


P3/Pi  =  (Us/U2)n 


10  log (P3/Prtf)  -  10  \og(Pj/ P ref)  =  n  10  \og[U3/U2) 


LP3  -  LPi  =  nlog(U3/U2) 


n  =  (LPi-LPi)/(\0\og(U3/U2)) 


where  the  subscripts  3  and  2  correspond  to  30  and  20  m/s.  Similar  equations 
hold  for  UtfU2,  U6/U2 ,  U3/U3  and  I/5/I/4.  Table  6  lists  the  calculated 

value  for  n  for  each  sample  tested.  The  values  for  n  varied  from  3.3  to  6.8  with 
a  mean  value  of  4.9  and  a  standard  deviation  of  0.3.  Excluding  the  extremely 
low  values  of  n  recorded  for  03  at  the  higher  free  stream  velocities,  the  mean 
value  of  n  was  5.0  with  a  standard  deviation  of  0.2. 
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The  speed/power  relationship  for  orifices  0 1,  02  and  louvers  L 1  and  L 2 
is  plotted  in  Figure  28.  Also  shown  is  the  least  squares  fit  to  the  data  points 
satisfying  the  equation: 


[LPi  -  LPj)  =  nlO  log {UJUj)  (5.2) 


From  the  fit,  the  best  value  for  n  is  4.97. 


Orifice/Louver 

U3/Ui 

UJUi 

Us/U 2 

u</u3 

tfs /U3 

Ub/U< 

0-1 

5.7 

5.4 

5.0 

5.2 

4.8 

4.5 

0-2 

6.0 

5.2 

4.2 

5.1 

4.4 

4.7 

0-3 

5.3 

4.5 

3.3 

4.2 

3.4 

3.5 

L-l 

4.6 

4.9 

5.4 

5.0 

5.4 

5.5 

L-2 

4.9 

5.2 

5.5 

5.1 

5.3 

5.0 

0-2A 

5.6 

5.1 

4.3 

5.1 

4.7 

5.2 

L-1A 

5.9 

5.1 

3.9 

4.9 

4.2 

4.5 

L-2A 

5.5 

4.8 

3.8 

5.3 

5.1 

6.8 

average 

5.4 

5.0 

4.5 

5.0 

4.7 

5.0 

std.  dev. 

0.5 

0.3 

0.6 

0.3 

0.6 

0.9 

overall  average  4.9 

overall 

standard 

deviation 

0.3 

Table  6:  Speed/Power  relationship  [P  =  CUn),  values  for  the  speed  exponent,  n 


Thus  the  concluding  dependence  is: 

Total  radiated  power  increases  with  the  5th  power  of  the  free  stream  velocity. 

5.8  Directivity  of  orifice  Ol  and  louver  L2A 


Directivity  measurements  from  orifice  01  and  louver  L2A  were  taken  at 
U  =  25  m/s.  Figure  24  details  the  coordinate  system  used  for  the  measurements. 
The  measurements  indicated  the  source  was  nearly  omni-directional.  Power 
levels  in  the  XZ  (azimuth)  plane  deviated  by  0.9  dB  from  the  maximum  level. 
In  the  YZ  (polar)  plane,  power  levels  deviated  by  2.6  dB  from  the  maximum 
level.  Off  the  major  axis,  power  levels  deviate  by  3.6  dB  from  the  maximum  level 
directly  above  the  orifice  source.  Table  7  lists  the  measured  sound  intensity  in 
dB  re  1 E  -  12  watt/ m2  at  a  radial  distance  of  r  =  27  inches  and  the  various 
polar  and  azimuth  angles.  Figure  35  is  the  azimuth  directivity  plot  for  orifice  01 
and  louver  L2A  and  clearly  supports  the  conclusion  of  a  nearly  omni-directional 
source. 
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Figure  33:  Speed/Power  relationship  for  four  samples  ( P  =  CUn )  with  least 
squares  fit 


Polar  angle 

Aiimuthal  angle 

Total  Intensity 

Orifice  0-1 

dBrelE  —  12  watt/m 2 

0 

90 

57.3 

0 

45 

56.9 

0 

15 

56.4 

0 

135 

57.1 

0 

165 

56.8 

0 

105 

57.2 

0 

75 

57.0 

-85 

90 

54.9 

-45 

90 

56.5 

0 

90 

57.5 

45 

90 

57.0 

85 

90 

55.9 

-45 

30 

54.7 

-45 

25 

53.9 

-45 

35 

55.0 

Louver  L-2A 

0 

160 

66.2 

0 

145 

66.1 

0 

130 

66.0 

0 

115 

66.0 

0 

90 

65.8 

0 

75 

65.7 

0 

60 

65.7 

0 

45 

65.8 

0 

30 

65.2 

0 

15 

65.2 

-90 

90 

63.4 

-75 

90 

63.9 

-60 

90 

65.0 

-45 

90 

65.4 

-30 

90 

65.9 

-15 

90 

66.1 

0 

90 

66.2 

15 

90 

66.4 

30 

90 

66.5 

45 

90 

66.6 

60 

90 

66.4 

75 

90 

66.3 

90 

90 

65.6 

-45 

45 

65.3 

-45 

30 

64.2 

-45 

90 

65.3 

-45 

135 

66.0 

Table  7:  Directivity  measurement  for  orifice  0-1  and  louver  L2A 
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Figure  34:  Coordinate  system  for  directivity  measurements.  Reference  is  the 
aperture  center.  Polar  angle,  6,  measured  from  Z  axis,  azimuthal  angle,  <f>,  mea¬ 
sured  from  —x  axis.  Measurements  described  by  spherical  coordinates,  r,6,4>. 

5.9  Edge  Angle  Effect  on  Radiated  Sound  Power 

5.9.1  Introduction 

Earlier  it  was  shown  that  indeed  the  orifice  source  provided  both  the  spectral 
character  and  the  predominant  amount  of  the  total  power.  The  strength  of  this 
source  is  due  in  part  to  the  eddy  interaction  with  the  downstream  edge  of  the 
orifice.  To  determine  the  relative  importance  of  this  edge  interaction  compared 
to  the  other  significant  parameters  identified  in  Section  3,  three  additional  ori¬ 
fice/louver  samples  were  tested.  They  were  orifice  0-2A  and  louvers  L-1A  and 
L-2A.  The  details  of  these  samples  are  described  in  Section  4.  The  only  differ¬ 
ence  between  these  samples  and  the  previously  tested  samples  was  the  leading 
and  trailing  edges  were  slanted  away  from  the  flow  by  40  degrees.  This  provided 
both  a  blunt  orifice  trailing  edge  where  the  eddies  strike  and  a  sharpened  orifice 
leading  edge  (essentially  reducing  the  upstream  aperture  to  1/32",  the  radius  of 
the  edge). 

5.9.2  Angled  edge  orifice/louver  power  spectra 

The  measured  sound  power  spectra  for  the  above  samples  ar5  shown  in  Fig¬ 
ures  31,  32,  and  33.  The  total  power  radiated  by  each  sample  at  the  various 
freestream  velocities  is  given  in  Table  5. 
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Figure  35:  Orifice  01  directivity  in  the  XZ  (azimuth)  plane  and 
directivity  in  the  YZ  (polar)  plane 
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Figure  39:  Edge  angle  effect  on  sound  power  spectral  level,  Ll  versus  L1A  at 
U  =  20  m/s 

Several  observations  are  worth  noting: 

•  the  effect  of  angling  the  edges  is  more  predominant  for  the  louvers  than 
the  single  orifice. 

•  the  effect  on  louver  L-1A  is  a  14-16  dB  reduction  in  total  radiated  sound 
power.  The  large  resonant  peaks  which  occurred  in  the  spectra  for  L-l  are 
not  present  in  the  spectra  for  L-l  A.  Figures  39  and  40  compare  the  two 
louvers  at  U  =  20  and  50  m/s.  Not  only  are  the  resonant  peaks  eliminated 
by  angling  the  edges  but  also  the  radiated  power  is  reduced  by  5-7  dB  in 
the  higher  frequency  bands. 

•  Figure  32  compares  the  power  spectral  levels  of  L-2  and  L-2A.  The  result 
for  this  louver  is  a  8.6  to  10.4  dB  reduc  tion  in  sound  power  equally  divided 
across  the  frequency  bands. 
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Figure  41:  Edge  angle  effect  on  sound  power  spectral  level,  L2  versus  L2A  at 
U  —  50  m/s. 

•  for  the  single  orifice,  0-2A,  the  effect  of  angling  the  edges  resulted  in  only 
a  0.7  to  1.4  dB  reduction  in  sound  power.  Figures  33  and  34  compare  the 
power  spectra  for  orifice  0-2  and  0-2A  at  U  =  20  and  50  m/s.  Here  again 
the  effect  of  angling  the  edges  resulted  in  a  reduction  in  sound  power  in 
the  higher  frequency  bands  but  in  this  case,  essentially  no  change  in  the 
lower  frequency  bands. 

5.9.3  Combined  effect  of  louvering  and  angling  edges 

Taken  together,  the  effect  of  louvering  and  angling  edges  results  in  a  reduc¬ 
tion  of  3.4  to  7.5  dB  in  radiated  sound  compared  to  the  single  angled  orifice 
(both  of  equal  total  open  area).  This  reduction  most  probably  is  caused  by  two 
factors: 
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Figure  44:  Effect  of  Louvering  on  sound  power  spectra  level,  L-2A  versus  0-2A 
at  U  =  30  m/s 

•  the  louvers  shed  eddies  which  after  convecting  downstream  to  the  trailing 
edge  are  the  same  order  in  size  as  the  edge  thickness.  This  results  in  a 
more  beneficial  interaction  when  the  eddies  impact  the  angled  edge.  The 
orifice,  however,  with  a  four  times  larger  aperture,  sheds  eddies  which  are 
much  larger  in  size  upon  impact  with  the  trailing  edge  and  are  therefore 
less  affected  by  the  thin  angled  edge. 

•  louvering,  by  reducing  the  transverse  dimension,  6,  interferes  with  the 
longitudinal  correlation  of  the  shed  eddies.  This  interference  will  be  more 
predominant  in  the  lower  frequencies  where  the  longitudinal  correlation 
length  is  larger  than  the  aperture  transverse  dimension,  b.  As  shown  in 
Figures  35,  36,  37,  and  38,  the  reduction  in  sound  power  is  predominantly 
in  the  lower  frequency  bands.  There  is,  however,  more  radiated  power  in 
the  angled  louver  in  the  higher  bands  by  approximately  1-3  dB. 
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Thus  one  can  conclude  that  effective  louvering  does  reduce  the  total  radiated 
power  by  reducing  single  mode  resonance  and  by  inter  fering  with  the  lower 
frequency,  longitudinal  eddy  correlation.  This  results  in  effectively  broadening 
the  power  spectrum. 

5.9.4  Louver  geometry  effect  on  sound  power 

As  with  the  straight  walled  louvers,  the  multiple  circular  aperture  louver 
was  quieter  than  the  multiple  rectangular  aperture  louver  by  1.4  to  3.6  dB 
in  total  radiated  power.  As  shown  in  Figures  39  and  40,  louver  L-2A  had 
approximately  the  same  spectral  shape  as  louver  L-l  A  but  with  considerably  less 
power  in  the  low  frequency  bands.  Again  this  reduction  can  be  attributed  to  the 
additional  interference  of  the  lateral  correlation  due  to  the  reduced  dimension  in 
the  aperture’s  height.  To  test  this  conclusion,  one  can  calculate  the  frequency 
where  the  lateral  correlation  length,  t ,  equals  the  aperture  diameter,  d  for  louver 
L-2A  and  at  freestream  velocity  of  50  m/s.  From  Figure  8,  at  B(Kr3)  —  0.4, 
the  value  of  Kr3  equals  1.2.  Setting  Kr3  equal  to  0.01  m  (27/64")  and  U  — 
50  m/s  gives  a  limiting  frequency  of  450  Hz.  Thus  one  would  predict  diminishing 
reductions  in  sound  power  for  frequencies  above  450  Hz.  Reference  to  Figure  49 
shows  that,  indeed,  at  frequencies  above  450  Hz  the  effect  of  lateral  correlation 
interference  on  reducing  sound  power  is  diminishing. 

5.10  Describing  the  Velocity  Field  Surrounding  the  Samples 

In  Section  3,  the  radiated  sound  power  was  shown  to  be  directly  related 
to  the  normal  component  of  the  turbulent  intensity.  Several  measurements  of 
the  turbulent  intensity  spectrum  were  taken  as  described  in  Section  4.4.  The 
assumption  was  made  that  the  turbulent  intensity  in  the  streamwise  direction 
was  six  times  that  in  the  normal  direction,  based  upon  measurements  of  corre¬ 
sponding  r.m.s.  velocity  component  values  in  the  most  active  region  of  turbulent 
boundary  layer  near  the  wall,  Klebanoff  [18].  Thus  8  dB  was  subtracted  from 
the  measured  intensity  spectral  levels  to  give  only  the  normal  component  of  the 
turbulent  fluctuations.  The  following  observations  are  made: 

•  the  turbulent  intensity  measured  2  tenths  of  an  inch  into  the  boundary 
layer  was  relatively  unaffected  by  the  probe  location  in  the  transverse 
direction  and  the  various  orifice  geometries.  Figure  50  is  the  measured 
normal  component  of  the  turbulent  intensity  for  louver  L-l,  with  the  probe 
at  the  center  of  the  upstream  aperture.  These  spectral  levels  were  used  in 
the  following  section  to  collapse  the  sound  power  spectra. 
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•  there  were  significant  reductions  in  turbulent  intensity  (approximately  3 
dB  )  between  louvers  L-l  and  L-1A  when  measured  at  the  leading  edge  of 
the  leading  aperture.  This  reduction  in  intensity  supports  the  measured 
reduction  in  total  sound  power  and  also  supports  the  conclusions  that  the 
sound  power  is  related  to  the  driving  aerodynamic  field. 

Despite  a  fairly  rigorous  attempt  to  maintain  the  block-house  air  tight,  a 
finite  and  measurable  pressure  differential  existed  between  the  blockhouse  and 
the  duct  static  pressure.  This  pressure  differential  resulted  in  a  small  ampli¬ 
tude,  mean  flow  entering  the  orifice/louver  from  the  closed  acoustic  space.  The 
magnitude  of  this  transverse  mean  flow  was  greatest  in  the  louver  samples  and 
measured  on  average,  fourteen  percent  of  the  freestream  velocity.  The  exact 
contribution  of  this  mean  flow  to  the  radiated  sound  power  in  not  known  but  it 
is  estimated  to  be  of  minimal  impact  since  the  turbulent  intensity  of  the  entering 
flow  measured  several  orders  of  magnitude  less  than  the  normal  component  of 
the  freestream  turbulence. 

5.11  Verifying  the  Scaling  Laws 

In  Section  3,  an  equivalent  source  model  was  developed  for  rectangular 
shaped  aperture  orifices  and  louvers.  To  ascertain  the  validity  of  this  model, 
the  sound  power  spectra  for  orifices  0-2,  0-3  and  the  louver  L-l  were  collapsed 
by  normalizing  the  spectral  levels  by  the  magnitude  of  the  equivalent  source 
model.  The  resulting  collapsed  spectra  are  shown  in  Figures  51,  52  and  53.  It 
is  evident  by  the  collapse  that  the  equivalent  source  model  adequately  describes 
the  controlling  parameters  for  orifice  and  louver  flow  noise.  The  prediction  of 
the  actual  magnitude  of  the  sound  power  is  fairly  accurate  for  the  samples  with 
the  smaller  apertures  but  is  too  high  by  about  7  dB  for  the  larger  aperture 
openings  (0-1  and  0-2).  The  probable  reasons  for  this  are  (1)  overestimating 
the  effective  aperture  width,  b,  and  (2)  the  estimation  for  the  lateral  correlation 
lengths.  As  mentioned  in  Section  3,  for  orifice  0-2,  the  longitudinal  correlation 
amplitude  function,  A(K(cj))  does  attenuate  to  a  value  of  approximately  0.6  for 
the  slower  speeds  indicating  that  some  eddy  decay  occurred  in  the  aperture.  For 
this  case,  one  could  assume  a  smaller  transverse  length,  6',  over  which  the  eddy 
decay  was  minimal  and  use  this  length  scale  for  evaluating  the  source  strength. 
Defining  6'  at  the  point  where  the  magnitude  of  the  amplitude  function  equaled 
0.95  results  in  an  effective  aperture  length  of  0.5  inches  versus  the  actual  dimen¬ 
sion  of  6  =  1.67  inches  (for  orifice  0-2  at  U  =  20  m/s  and  /  =  250  Hz).  This 
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Figure  51:  Scaling  law  verification  for  orifice  0-2 


results  in  a  reduction  of  4.2  dB  in  the  predicted  spectral  level  and  could  account 
for  some  of  the  error  in  predicting  the  larger  aperture  radiated  sound  power. 
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6.  CONCLUSIONS 


The  acoustic  source  created  by  turbulent  flow  over  orifices  and  louvers  was 
investigated  by  measuring  the  radiated  sound  intensity,  utilizing  an  intensity 
probe,  at  various  wind  speeds  in  a  semi-anechoic  wind  tunnel.  From  the  intensity 
spectra,  the  radiated  power  and  the  directivity  of  the  source  was  determined. 
The  goal  of  this  research  was  to  describe  the  source  and  identify  the  controlling 
parameters  which  affect  its  spectral  shape,  radiated  power  and  directivity.  The 
orifice  and  louvers  investigated  have  geometries  similar  to  those  used  in  industry 
in  both  air  and  water  environments.  The  effects  of  various  leading/trailing  edge 
geometries  were  also  investigated. 

The  source  of  noise  has  been  modelled  on  the  basis  of  theoretical  analyses 
by  Ffowcs  Williams,  Nelson  and  Corcos.  The  equivalent  source  of  noise  was 
modelled  as  two  equal  but  opposite  monopole  sources  located  in  the  upper  and 
lower  aperture  of  each  orifice.  The  strength  of  each  source  is  determined  by 
the  mass  flow  driven  through  the  apertures  by  the  normal  component  of  the 
turbulent  fluctuations.  The  resulting  radiated  sound  field,  both  into  a  closed 
acoustic  space  and  radiating  away  from  the  orifice,  is  dependent  on  the  acoustical 
properties  of  the  louver. 

To  ascertain  the  validity  of  the  source  model,  two  orifice  and  one  louver  power 
spectra  were  collapsed  by  normalizing  the  spectral  levels  by  the  magnitude  of  the 
equivalent  source  model.  It  is  evident  by  the  collapse  that  the  equivalent  source 
model  adequately  describes  the  controlling  parameters  for  orifice  and  louver  flow 
noise.  The  controlling  parameters  for  the  geometries  tested  were: 

1.  Strouhal  number  based  on  the  lateral  correlation  length,  f,  and  aperture 
transverse  dimension,  6 

2.  the  normal  component  of  the  Turbulent  Intensity, 

3.  the  fifth  power  of  the  freestream  velocity  and 

4.  the  shape  factor  (1  -  cos(/C6))  where  K  is  the  convective  wavenumber, 
K  =  u/Uc,  and  Uc  is  the  eddy  convection  velocity. 

As  a  result  of  this  research,  several  important  conclusions  are  made: 

•  sound  intensity  measurements,  using  an  intensity  probe,  can  successfully 
describe  an  acoustic  source,  even  though  measurements  are  taken  within 
the  source’s  near  field. 


•  the  errors  associated  with  measuring  intensity,  using  an  intensity  probe, 
can  be  as  small  as  the  inherent  statistical  deviation  of  the  random  fluctu¬ 
ations. 

•  the  equivalent  source  model  developed  to  describe  the  behavior  of  ori¬ 
fice  and  louver  flow  noise  successfully  describes  the  controlling  parameters 
and  provides  valuable  insight  for  predicting  the  actual  magnitudes  of  the 
radiated  sound  spectrum. 

•  the  1/3  octave  band  containing  the  most  acoustic  power  from  orifice  and 
louver  flow  noise  can  be  predicted  from  the  non-dimensional  Strouhal  num¬ 
ber  based  on  the  aperture’s  transverse  dimension  and  the  non-dimensional 
length  scale, 6/6  where  6  is  the  turbulent  boundary  layer  thickness. 

•  for  aerodynamic  flows,  the  orifice/louver  source  provides  both  the  domi¬ 
nant  amount  of  total  sound  power  and  the  spectral  character  of  the  ra¬ 
diated  noise.  Contributions  due  to  structural  radiation  from  the  ducting 
and  ducting  “open  window”  effects  are  not  significant. 

•  the  effect  of  louvering,  simply  by  dividing  the  total  open  area  into  more 
numerous,  smaller  areas  does  not  by  itself  reduce  the  total  radiated  flow 
noise  and  will  most  probably  result  in  increased  radiated  sound  power  in 
the  higher  frequency  bands. 

•  the  effect  of  reducing  the  aperture  throat  thickness  is  to  broaden  the 
spectral  shape  of  the  radiated  power  spectrum  while  maintaining  approx¬ 
imately  the  same  total  radiated  sound  power. 

•  the  effect  of  reducing  the  aperture’s  transverse  dimension  can  be  ade¬ 
quately  described  by  the  controlling  shape  factor  in  the  equivalent  source 
model,  (l  -  cos(Kb)). 

•  directivity  measurements  show  a  nearly  omni-directional  source  in  the  on- 
axis  planes  with  off-axis  power  levels  reduced  by  2-3  dB. 

•  taken  together,  the  effect  of  louvering  and  angling  the  leading  and  trailing 
edges  resulted  in  reductions  to  sound  power  of  3.4  -7.5  dB  compared  to 
the  equivalent  area,  single  aperture,  angled  orifice. 

•  the  measured  reductions  in  sound  power  due  to  effective  louvering  can  be 
explained  by  (1)  reductions  in  the  normal  component  of  the  turbulent  in¬ 
tensity  due  to  the  eddies  impacting  a  blunt  trailing  edge  and  (2)  destructive 
interference  with  the  eddy  longitudinal  and  lateral  correlations. 
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the  louvers  with  numerous  circular  apertures  and  angled  edges  radiated  j 

the  least  acoustic  power  with  a  very  broad,  flat  power  spectrum. 
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APPENDIX  1: 

I  MATCHED  ASYMPTOTIC  EXPANSION  SOLUTION 

TO  A  RADIATING  ORIFICE  IN  AN 
INFINITE  RIGID  WALL 

I  The  two  dimensional  problem  to  be  solved  is  the  far  field  radiated  pressure 

at  I*,  created  by  a  driving  aerodynamic  field  at  yt,  within  a  finite  size  aperture 
in  an  infinite,  rigid  wall.  Using  the  principle  of  reciprocity,  the  above  field  can 
be  determined  by  first  solving  for  the  sound  field  at  yt  created  by  a  distant 
monopole  at  it-  Then,  the  desired  solution  is  obtained  by  simply  replacing  the 

S  position  vectors  of  the  far  field  monopole  and  the  aperture  source.  The  solution 

to  this  problem  will  be  analyzed  using  the  principles  of  Matched  Asymptotic 
Expansions  (MAE). 

The  coordinate  system  used  in  the  analysis  is  shown  below. 


r. 


with  boundary  conditions,  dp/dy j  =  0  for  J  yt  |>  b  and  y*  =  0. 

The  governing  equations  for  a  plane  wave  of  single  frequency,  u,  at  normal 
incidence  are: 


V2p  +  fcJp  =  0 


(A  1-1) 


p  =  Ae'[ki'-Ut) 


(.41.2) 
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where  k  is  the  acoustic  wavenumber,  k  =  ui/co.  Also  define  the  non-dimensional 
terms,  r  =  u/t,y i  =  ky 2,y2  =  ky2,p  =  p/A  and  e  —  kb  for  the  parameters  in  the 
outer  problem. 

Substituting  the  above  outer  parameters  into  equations  Al.l  and  A1.2  gives 
the  governing  equations  and  boundary  condition  for  the  outer  problem 


V’p  +  p  =  0 


(A1.3) 


p  = 


(A1.4) 


To  lowest  order,  e  — ►  0,  and  the  orifice  closes.  Thus  the  outer  problem 
becomes 


P  ~  ^iPi  +  62p2  +  •  •  •  as  e  — *•  0 


(/tl.5) 


with  =  1.  From  the  boundary  condition,  we  can  determine  pi  to  be 


Pi  =  P.  +  Pr  =  e‘w’  +  t  V3 ,  y-i  <  0 


(A1.6) 


suppressing  the  e~'r  factor,  where  p,,pr  are  the  incidence  and  reflected  waves, 
respectively.  Note  that  the  boundary  condition  is  satisfied 

=  o  for  y2  =  0.  (A1.7) 

dy2 


Choose  also  the  condition  that  pi  =  0  for  y2  >  0.  To  open  the  orifice,  choose 
the  inner  variables  yi  =  j h/t  and  y2  =  y2 /t. 

In  inner  variables, 
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P  =  exp[i(ey2  -  f)]  +  exp[-i{ty2  -  f)] 

=  exp(-if){exp((yt)  +  exp{~eyt)} 

fly*  e2S2 

~  ezp(-*r){l  +  ty2 - - - +  1  -  cya - ^+  --} 

£*t/* 

=  2e  ,r{l - — *  +  0(£4)}  as  c -*  0,y2  fixed  (y2  -»  0).  (>11.8) 

h 

The  governing  equation  with  the  inner  parameters  substituted  into  equations 
A  1.3  becomes 

V2p  +  £sp  =  0  (A1.9) 


with  boundary  conditions  dp/dy2  =  0  for  y2  =  0  and  |  |>  1.  Let  p  be  defined 

by 

p  ~  Aq(£)pi  +  Ai(f)p2  +  ■  •  •  (A1.10) 


for  the  inner  expansion.  The  resulting  inner  problem  coordinate  system  is 


y% 


Since  V2pj  =  0,  solve  the  inner  problem  by  conformal  mapping,  pi  ~ 
ARe(w).  Using  a  solution  of  the  form,  iz  =  coshu;,  (suppressing  the  t~'T  factor) 
we  can  determine 

11-3 

w  ~  tn2z  +  itt/2  H - - - 7  +  •  •  •  as  I  z  I— ►  oo  (-41.11) 

'  2  •  2z2  2  •  2  ■  4z4  11  v  7 
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where 


Since 


z  =  fe'#. 


z  =  z/e, 


u;  ~  /n2z  —  £ne  +  *?r/2  +  0(e:)  for  e  — ♦  0,  2  fixed(z  — ►  00)  (,41.12) 


Therefore,  utilizing  equation  A. 11,  we  can  determine  A0(e)pi  to  be 
A0(e)pi  ~  AA0{£n2f  -  tne  +  iir/2  +  0(e1)}. 


(A1.13) 


Now  matching  the  leading  terms  in  the  outer  and  inner  solutions,  equations 
A1.8  and  Al.13,  leads  to  A0  =  +1  /tnt  and  A  =  -2.  Then  <52(e)  =  l/lnt  in 
the  outer  expansion.  p2  satisfies  the  Helmholtz  equation  with  k  =  1  and  must 
behave  like  tnr  for  r  — >  0,  t  — >  0  and  f  fixed.  Clearly 


P2  =  CH{0'\r) 


(A1.14) 


satisfies  this  requirement  and  also  satisfies  the  rigid  wall  condition.  As  r  — ►  0, 
H^(r)  ~  —  tnr  +  1  +  —  (7  -  tn2)  +  •  •  • 

7T  7T 


where  7  is  defined  as  Euler’s  constant.  Thus  to  match  with  the  inner  solution, 
A  =  —  2  =  C  •  2»/jt,  C  must  be  equal  to  iir.  Therefore  the  final  outer  expansion 
becomes 
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\  lm  I 


(41.15) 


where  the  first  term  in  brackets  is  the  incident  wave,  the  second  term  is  the 
reflected  wave  and  the  lowest  order  orifice  effect  is  the  last  term  which  describes 
a  monopole  type  radiation.  Substituting  for  the  non-dimensional  terms  and 
applying  the  principle  of  reciprocity,  one  obtains  the  radiated  sound  from  a  2 
dimensional,  single  slot,  line  source 


P 


A 


t 


(41.16) 


As  r  — *  oo ,ffn^(f)  —*  so  the  far  field  radiated  sound  for  the 

2D,  single  slot,  line  source  becomes 

Pr-oo  -  A  {c4il  +  e"kil  +  t7r^/2/jr?e’(r'”#/4)|  e-,w<.  (41.17) 


Equation  A1.17  provides  the  result  that  the  far  field  radiated  pressure  for  a 
2D,  single  slot  source  to  lowest  order  is  composed  of: 

1.  the  incident  and  rigid  wall  reflection 

2.  2D  line  source  having  a  1/r  dependence  (cylindrical  spreading). 

Thus  one  can  infer  that  it  is  plausible  for  the  3  dimensional,  single  orifice 
source  to  have  a  reflection  coefficient,  |  R  j*  as  1  for  values  of  6/A  <c  1. 
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APPENDIX  2: 

MEASURING  TURBULENT  INTENSITY 

In  this  appendix,  the  transfer  function  to  relate  the  mean  and  fluctuating 
component  of  the  hot  wire  anemometer  output  voltage  to  the  mean  and  fluctu¬ 
ating  component  of  the  particle  velocity  is  developed. 

The  basic  relationship  between  output  voltage,  E,  and  particle  velocity,  U, 
is  defined  by  King’s  Law 

E  =  A  +  BUn  (A2.1) 


Now  define  the  instantaneous  voltage  and  particle  velocity  by 

E  =  E0  +  t'  and  |  e'  |<|  E0  \  (A2.2) 


U  —  i/0  +  u'  and  )  u'  |«:|  U0  j  (A2.3) 


where  E0  is  the  mean  of  E,  U0  is  the  mean  of  U,  e'  and  u'  are  the  fluctuating 
components  whose  expectations  are  zero.  Substituting  equations  A2.2  and  A2.3 
into  equation  A2.1,  expanding  and  dropping  the  higher  order  terms  of  u'  and  e' 
gives 

El  +  2  E0c'  =  A  +  BU£{1  +  [nu'/U0})  {A2A) 


Taking  the  expectation  of  equation  A2.4  gives 

El  =  A  +  BU£.  (A2.5) 
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Subtracting  equation  A2.5  from  equation  A2.4  gives  the  desired  relationship 
between  the  fluctuating  components  of  the  anemometer  output  voltage,  e',  and 
the  particle  velocity,  u', 

u'  =  {[2  E0Ul0-n\!  Bn)t' .  (A2.6) 


Thus,  defining  the  term  in  parentheses  as  H ,  the  resulting  relationship  be¬ 
tween  the  spectral  densities  of  the  fluctuating  velocity,  Su(/),  and  the  output 
voltage,  SE(f),  is 


APPENDIX  3: 

MATHEMATICAL  STEPS  IN  DETERMINING 

Sq(yk»w) 

From  equation  3.7,  the  cross  spectral  density  is  assumed  to  be  of  the  form 
Sw(f*,u>)  =  Sv{yk^)A{Kei)B(Kei)  •  ezp(iFfci).  (A3.1) 


As  described  in  Section  3.2,  the  longitudinal  correlation  amplitude  function, 
A(if«i),  is  identically  equal  to  one  for  all  samples  except  0-1  and  0-2  at  the 
lower  freestream  velocities,  where  at  a  minimum  it  equals  0.6.  So  assuming  no 
decay  in  the  eddy  structure  in  the  yi  direction  (i.e.,  equals  one)  and 

S„(y*,w)  to  be  uniform  over  the  aperture  surface,  the  source  strength  spectral 
density  becomes 

SQ{yK,u),=  p2u2Sv(yk,u)  f  f  5(/fcs)cip(:'A'f1)0(e*)dA(f/t)  (A3. 2) 

J  oo  J  oo 


where  ej,e3  are  the  position  coordinates  as  shown  below. 
Substituting  for  fl(e4)  we  have 


Sq{yk,u)  =  p2w2S„( y*,w)  f  f  B(/Ce3)e'*“(L- 

J  oo  J  OO 


f3  \)[b-  |  fi  I)  •  dA{ €k). 

(A.3.2o) 


Since  the  lateral  correlation  l  is  typically  one-tenth  or  less  the  length  of  the 
aperture,  the  integral  in  the  e3  direction  can  be  evaluated,  independent  of  the 
actual  limits  of  the  integrand,  as 

/  B(Jfe,)(L-|e,|)*,  =  W-0(O  (A3.3) 

J  OO 


where 
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Evaluating  the  integral  over  the  ti  directions  gives: 

eiktl  (6—  |  €,  |)dci  =  ^(1  -cos(tf6)).  (/13.4) 


Therefore  the  final  result,  neglecting  the  l1  term  is: 

S$(y*,w)  =  plwlSu(y*,w)L£(-^)(l  -  cos(Kfe)).  (43.5) 


The  remaining  term  needed  to  fully  describe  the  magnitude  of  the  source 
strength  spectral  density  is  t.  The  lateral  correlation  length,  l,  is  determined  at 
the  point  where  the  amplitude  function,  B[Kts)  equals  0.4  and  the  frequency 
equals  the  center  frequency,  fc ,  as  defined  in  Section  5.3. 


